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Abstract. Metamagnetism in layered ruthenates has been interpreted as a novel kind of quantum critical 
behavior. In an external magnetic field, Ca2-a;Sri,Ru04 undergoes a metamagnetic transition accompanied 
by a pronounced magnetostriction eff'ect. In this paper we present a mean-field study for a microscopic 
model that naturally reproduces the key features of this system. The phase diagram calculated is equivalent 
to the experimental T-x phase diagram. The presented model also gives a good basis to discuss the critical 
metamagnetic behavior measured in the system. 
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! 1 Introduction 

Ca2_xSrxRu04 has attracted interest for its complex and 
puzzling phase diagram including metallic as well as Mott- 
insulating magnetic phases pL^ which depend in a subtle 
' way on structural properties. Recent magnetostriction ex- 
, periments |2|3| in connection with the metamagnetic tran- 
■ sition (MMT) for a: « 0.2 underline the strong coupling of 
\ the electronic properties to the lattice. This may be taken 
. as a hint for the relevance of localized electronic orbital 
and spin degrees of freedom as can be found in a num- 
, ber of transition metal oxides. In particular, the mutual 
• influence of spin and orbital degrees of freedom plays an 
\ important role in the behavior of manganites, ruthenates 
I or titanates Hj. 

The result of the apparent interplay between orbital 
, and magnetic correlations in Ca2_2:Sr2:Ru04can be seen 
in the complex T-x phase diagram shown in Fig. [1] Both 
Sr^+ and Ca^"'' are isovalent ions so that the substitu- 
tion of one by the other does not change the number 
of conductance electrons. Sr2Ru04 is a good metal and 
the a priori expected change due to the substitution of 
Sr by Ca should be increasing metallicity, because the 
doping with a smaller ion (Ca) would imply a widening 
of the band. This is not the case because the smaller 
Ca-ion induces lattice distortions which alter the over- 
lap integrals and the crystal fields of the relevant elec- 
tronic orbitals. Actually, Ca2Ru04 behaves as a Mott in- 
sulator, and the evolution between these end-members 
builds a very rich phase diagram where different struc- 
tural and magnetic phases appear (Fig[T]). For x = 2, cor- 
responding to Sr2Ru04, the system has tetragonal sym- 
metry with the RuOe-octahedra slightly elongated along 



the c-axis. Ca substitution initially induces the rotation of 
RuOg octahedra around the c-axis in order to accommo- 
date the smaller ions. The system behaves still as a param- 
agnetic metal with tetragonal symmetry. Further doping 
with Ca reduces the conductivity and the susceptibility 
increases when approaching x 0.5, becoming a Curie- 
like susceptibility [S]. At x = 0.5 there is a structural 
phase transition, where the crystallographic structure of 
the Ca2-2:Sr2,Ru04 series changes from tetragonal to or- 
thorhombic through a second-order phase transition. For 
X < 0.5 there is, besides the c-axis rotation, a tilting of 
the RuOg octahedra leading to a reduction of the symme- 
try and a reduction of the c-axis lattice constant. More- 
over, these distortions are responsible for a narrowing of 
the conduction bands which in turn enhances the correla- 
tion effects. In the region 0.2 < x < 0.5 the experiments 
suggest the appearance of a low-temperature antiferro- 
magnetic (AFM) order. Furthermore, the application of 
a magnetic field leads to a metamagnetic transition which 
influences the tilting of the RuOg-octahedra 2 . Within 
the region x < 0.2 the system is a Mott-insulator with 
true long-range AFM order and a total spin of 5 = 1 . 

In general, a MMT can occur for a material which, un- 
der the application of an external fleld, undergoes a first 
order transition or is close to the critical endpoint of such 
a transition to a phase with strong ferromagnetic correla- 
tions. This is experimentally observed by a very rapid in- 
crease of the magnetization over a narrow range of applied 
magnetic field. The field dependence of the magnetization 
and magnetoresistance of Ca2-a;Sra;Ru04 has been ana- 
lyzed in Ref. [5]. There, a MMT to a highly polarized state 
with a local moment of 5 = 1/2 is found. These measure- 
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Fig. 1. Sketch of the temperature-doping phase diagram of 
Ca2-a;Sr2;Ru04. Region I corresponds to a Mott insulator with 
long range antiferromagnetic order. Region II is characterized 
by a metallic behavior with orthorhombic lattice symmetry and 
antiferromagnetic correlations at low temperatures. Region III 
corresponds to a paramagnetic metal with tetragonal symme- 
try in the lattice. For x 0.5 there are strong ferromagnetic 
correlations at low temperatures. 

ments are interestingly supplemented by the magnetostric- 
tion experiments published in Ref. 2 3 where the change 
of the lattice constants as a function of the magnetic field 
for different temperatures is shown. Namely, these results 
demonstrate that crossing the MMT leads to an elonga- 
tion of the c axis accompanied by a shrinking along both 
in-plane directions. Apparently, the application of a high 
magnetic field at low temperatures can reverse the struc- 
tural distortion that occurs for 0.2 < x < 0.5 upon cool- 
ing in zero field. Furthermore, apart from different energy 
scales, the qualitative effects associated with the MMT 
are independent of the magnetic field direction. 

All the members of the Ca2-2;Srj;Ru04 family have 
4 electrons in the t2g orbitals of the Ru 4(i shell. What 
is different is the occupation of these orbitals in the two 
end- members of the phase diagram: while for a; = there 
is an average occupation of two electrons in the dxy orbital 
and the other two in the dyz and dzx orbitals, for x = 2 
there is a fractional occupation of 4/3 in the dxy band 
and 8/3 in the dj^^-d^jj-bands. LDA calculations for this 
concentration give three Fermi surface sheets, one with 
essentially xy and two with mixed {xz,yz} character [7]. 
The first one is usually called 7 band while the others are 
labeled by a and /? bands. For intermediate values of x the 
orbital occupation is still a matter of debate [8|9|10|ll|12j . 

In this work we focus on the region II and the Ca- 
rich part of region III in the schematic T — x phase dia- 
gram. For the microscopic description we follow the sce- 
nario of an orbital-selective Mott insulator. This scenario 
was put forward by Anisimov et al. [H] in order to ex- 
plain the unexpected effective magnetic moment close to 
a 5 = 1/2 spin for 0.2 < a: < 1.5 [6]. Assuming that 
the orbital occupation in this region is (na^p,n~^) « (3, 1), 
they proposed that the electrons in the {a, /?} bands un- 
dergo an orbital-selective Mott transition (OSMT) while 
the 7 band remains metallic. In this scenario the exper- 
imental observation of the 1/2 effective spin is assigned 



to the localized hole in the {a, /?} bands. Angular magne- 
toresistance oscillations measurements jl3| indeed show a 
strong dependence of the Fermi surface on the Ca con- 
centration which is consistent with the scenario of co- 
existing itinerant and localized d-electronic states. Fur- 
thermore, from the theoretical point of view, there is by 
now a consensus that an OSMT can in principle occur in 
multi-band Hubbard models under rather general condi- 
tions [14115116117118119120121122] . However, it is stiU un- 
clear to which extend this concept is applicable in the 
Ca2--2;Sra:Ru04 system 123 24 . Despite of these uncer- 
tainties, there is little doubt that the localized degrees of 
freedom play an important role for the understanding of 
the puzzling physics of this material, and in the following 
we will assume that the concept of the OSMT is valid as a 
lowest order picture for Ca concentrations corresponding 
to region II and at the boundary of region HI of the T — x 
phase diagram. 

Using a mean-field description we focus on the inter- 
play between structural distortion and magnetic and or- 
bital ordering. We find a theoretical phase diagram which 
can be related to the experimental one. Furthermore, our 
calculations qualitatively reproduce the metamagnetic tran- 
sition accompanied by a structural transition observed in 
the system. The paper is organized as follows: in Sec. [2] 
we introduce the microscopic model. A mean-field analy- 
sis is performed in Sec. [3] The results for zero and finite 
magnetic field are given in Sec. [H In Sec. [5] we relate our 
results to the experimental measurements. We summarize 
the main conclusions of this work in Sec. [S] 

2 The model 

Following the scenario of the OSMT we assume one hole 
in the {a, /?} bands and focus on the localized orbital and 
spin degrees of freedom |8l25j . Neglecting for the our dis- 
cussion the itinerant 7 band it is natural to consider a 
two-dimensional extended Hubbard model of the form 

i,a,s 

+U ^ ^ ni^-^rii^i + ^' X! "•^zx'^i.yz 

i L' i 
'^Jh ^\,yz,s^\,zx,s''^i,^^,sCi^yz,s' (1) 

where c\ ^ {ci.v,s) creates (annihilates) an electron on site 
i with orbital index v (= yz, zx) and spin s {ui^^^s ~ 

4,iy,sCi,i/,s, = ?^^,l/,T a = {ax.ay) = (1,0) or 

(0, 1) basis lattice vector). With this Hamiltonian we re- 
strict ourselves to nearest-neighbor hopping and on-site 
interaction for the intra- and inter-orbital Coulomb re- 
pulsion, U and [/', respectively, and the Hund's rule cou- 
pling Jh- The hopping terms considered in this model 
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come from the 7r-hybridization between the Ru-d and O- 
p-orbitals and lead to the formation of two independent 
quasi-one-dimensional bands: the band associated to the 
(ij,2-orbital disperses only in the y-direction while the band 
associated to the dzx-OThital disperses in the x-direction. 

In the strongly interacting limit it was proposed that 
the a-/3-bands absorb 3 of the four electrons available per 
site and form a Mott-insulating state with localized de- 
grees of freedom, spin 1/2 and orbital [8j. The local orbital 
degree of freedom can be represented as an isospin config- 
uration |-t-) and |— ) corresponding to the singly occupied 
dzx and dyz orbitals, respectively. The isospin operators 
therefore may be defined as: 

r|±)^±i|±>, /+!-) = |+), /-|+) = |->. 

Taking into account the additional spin 1/2 degree of free- 
dom (I 1) and I I)) leads to four possible configurations at 
each site, represented by the states 

{|T+), IT-), U+>, U->}. 

Within second order perturbation in </[/ it is possible to 
derive from Ti.a,f3 an effective model describing the inter- 
action between the localized degrees of freedom. One finds 
the following Kugel-Khomskii-type model [8]: 

+ [C{IU. + V.){i:+V'j + D]} (2) 

where J — At'^/U. We have imposed the approximatively 
valid relation U = U' + 2Jh and have assumed that a = 
U'/U > 1/3. The parameters A, B,C, D,r]a_ and 77^ are 
functions of a alone and have been given elsewhere |8|25j . 
The energy scale JC > of the isospin coupling is the 
largest in the present Hamiltonian. Therefore, in a mean- 
field approximation, one expects antiferro-orbital (AFO) 
order below a critical temperature Tafo ~ JC (C > 0). 
On the other hand, the value of the spin-spin interac- 
tion depends on the orbital order and lies between Ji = 
J[A{r]l - 1/4) + B] <0 and J2 = J[Ar]l + B]. Thus, in 
the presence of AFO order the spin will align ferromag- 
netically (FM) below a critical temperature TpM ^ —Ji- 
If, however, AFO order is suppressed, as in the case of 
an orthorhombic distortion (see below) the spin-spin cou- 
pling is given by J2. We mention here that the sign of 
J2 depends on the value of a. In particular, J2 < for 
a < ac = 0.535 and consequently we expect FM order at 
low temperatures whereas for a > ac we have J2 > and 
antiferromagnetic (AFM) order sets in at sufficiently low 
temperatures. To be consistent with experiments we will 
choose throughout this article a value a — 0.75. 

The Sr substitution for Ca acts as an effective nega- 
tive pressure. In order to account for the orthorhombic 
distortion due to the tilting of the RuOg octahcdra, we 
introduce a new term in the Hamiltonian Tidist, defined 
as 

ndrst = lGN{e-eo)^+KeJ2l"^ (3) 



where G is the elastic constant, N is the number of Ru 
atoms and K is a coupling constant, e is a strain-field 
which accounts for an orthorombic distortion. Any or- 
thorombic distortion yields a uniform bias for the local 
orbital configuration which suppresses AFO ordering. This 
is modeled by the coupling of e to the orbital degrees of 
freedom. In other words, the orthorombic distortion intro- 
duces a transverse field which aims to align the isospins. 
The first term in Eq. ([3]) is a measure of the lattice elastic 
energy. In addition to the strain e driven by orbital corre- 
lations we assume a constant contribution Eq. We do not 
specify further the origin of this contribution but it might 
include effects of the 7 band or other, non-electronic, mech- 
anisms. For actual calculations we fix the value at Eq = 
0.1. We will discuss later to which extend we can relate 
the elasticity G in the theoretical model to the Sr concen- 
tration X in Ca2-xSra;Ru04. 

Finally, in order to study the metamagnetic transition, 
we introduce a coupling of the system to a magnetic field, 
by the inclusion of the term Timag, 

nr,^ag^-9iiBHY,Sf, (4) 

i 

where g is the electron gyromagnetic factor, /is is the 
Bohr magneton, fu-B = 27^ ^'^d H is the magnetic field 
strength. With all these ingredients, the full Hamiltonian 
can be written as 

H = Heff + Hdist + Ti-mag- (5) 

In the next section we treat this model in a mean-field 
approximation. 

3 Mean-field analysis. 

The mean- field decoupling for Ti^f / reproduces well some 
of the experimentally observed features of Ca2-xSr2;Ru04 in 
the X region where the band filling corresponds to the 
(n(c_^),n^) = (3,1) orbital occupation, as shown in Ref. 

Here we extend this analysis to the full Hamiltonian 
Eq. ([5]) and obtain the zero and finite magnetic field phase 
diagrams where the different competing orders of the sys- 
tem are represented. In the presence of a transverse mag- 
netic field, there appears an uniform component of the 
magnetization in the x-direction, 

(5f)=mo, (6) 

as well as a staggered component in the z-direction, 

(^^)-{^. ;[ \Ib. 

Here, we have made use of the bipartite structure of the 
Hamiltonian Eq. ([5]): A and B label the two sublattices. 
In addition, we introduce the staggered isospin component 
in the z-direction 
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The partition function of the system can be calculated 
by Z{P) = Tre-^'^, where f3 = l/fcsT. On the mean-field 
level, the bipartite nature of the lattice splits the system 
into two subsystems, so we can express the partition func- 
tion as 



N 



n 



(9) 



i£A,B 



where Eq denotes the energy density (per hole) corre- 
sponding to the term in the mean-field Hamiltonian that 
does not couple to any spin or isospin operator, 



Eo = J 



QAtliml ~ ml) - 2iAfjl + B){ml - ml) 



+2Ctl + 2C7j'^ +2D 



+ ^ie-eo)^ 



(10) 



N = Na + Nb is the total number of sites. Now we can de- 
fine the orbital Zf and magnetic Zf one-particle partition 
functions as l26l: 



(11) 

(12) 



where i?* and are the molecular field vectors that 
couple to the isospin and spin degrees of freedom. If we 
denote the total free energy of the system by F, the free 
energy per site (per hole) T ~ F/N = — \nZ//3N is 

1 

2P 



T 



^^0 - — (In + In Z|j + In Z\ + In Z%) (13) 



Eventually, the free energy of the system for the given 
mean fields is found to be 



^Atliml - ml) - 2(^77^ + B){ml - ml) 



-2Ctl 
1 



2Crj'^- 



2 cosh 



2 cosh 



2D 



jie-so)' 



(14) 



where 



A = i^V ^ ^Qj2 [Atsimf, - ml) + Cts] , 

/2 - {-gfiBH + AJmo [A{r,l - tl) -f S] toq}' 

+ 16J'ml[A{rjl~tl) + B]\ (15) 

The values of the mean fields are determined by the solu- 
tion of the self-consistency equations 



(16) 



where (j) = £,ts,mo,ms. These equations are solved nu- 
merically. Because 



in 



G 
K 



(e - £o) 



the ferro-orbital order parameter is directly related to the 
difference e — Eq. 



01 




Fig. 2. Sketch of the low-temperature competing orders of the 
system. On the left hand side it is represented a lattice with 
orthorhombic symmetry and ferro-orbital order plus antifer- 
romagnetism (FO & AFM). On the right hand side, a lattice 
with tetragonal symmetry with antiferro-orbital order plus fer- 
romagnetism (AFO & FM). 



4 Results 

The competing low temperature orders are schematically 
shown in Fig. [2] For small values of G, the free energy min- 
imization gives a ground state with a strong lattice distor- 
tion breaking the tetragonal symmetry, as shown on the 
left hand side of Fig. [2l In this case, the symmetry of the 
lattice is orthorhombic and ferro-orbital order (FO) co- 
exists with antiferromagnetic order (AFM). On the other 
hand, for large values of the elastic constant G, anti-ferro 
orbital (AFO) and ferromagnetic (FM) spin order are si- 
multaneously realized. This is sketched on the right hand 
side of Fig. [1 



4.1 Absence of magnetic field. 

The G-T phase diagram for = is shown in Fig.[3l The 
elastic constant G of the lattice controls the distortion. 
Apart from the high-temperature disordered phase (PO) 
we can distinguish two main regions in the phase diagram, 
as discussed below. 



4.1.1 Soft lattice 

Lowering the temperature from the disordered phase we 
find for a soft lattice (small G) a crossover to a ferro- 
orbital ordered (FO) state. This crossover takes place at 

a temperature in the range of ksTpo = jq indicated by 
the diffuse line in Fig. [3] The FO order is accompanied 
by a substantial orthorhombic distortion e ^ Eq which is 
driven by the coupling of the strain field to the orbital 
degrees of freedom as described by Eq. ([3]). In addition, 



below a critical temperature ksTAFM = J {B 
tiferromagnetic (AFM) order sets in. 



Mi 



(17) 4.1.2 Hard lattice 



For a harder lattice (large G), the gain in energy by po- 
larizing the orbital degrees of freedom is not sufficient to 
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Fig. 3. T — G phase diagram for H — involving antiferromag- 
netic (AFM), ferromagnetic (FM), ferroorbital (FO) and anti- 
ferroorbital (AFO) order. Dashed lines indicate second order 
phcise transitions while full lines represent first order transi- 
tions. TP is a tricritical point and BP is a bicritical point. The 
diffuse line between FO and paraorbital (PO) region indicates 
a crossover. 



drive a substantial orthorhombic distortion and e ~ Eq. In- 
stead, there is a staggered orbital order (AFO) which sets 
in at temperatures slightly below ksTAPO = JC. The 
AFO order drives a ferromagnetic (FM) ordered phase 
roughly below ksTpM — — J [B + A (^ijl^ - j)]. For smaller 
values of G, this transition is weakly first-order but changes 
its character at the tricritical point (TP) to second-order. 



4.1.3 Structural transition 

There is a structural transition between the soft and the 
hard lattice. Below Tafm this transition is of first-order 
and is characterized by a simultaneous discontinuity in all 
the order parameters. The first-order line splits into two 
second-order lines at a bicritical point (BP). For temper- 
atures above Tap m the transition is characterized by the 
onset of a staggered orbital and a gradual suppression of 
the ferro-orbital order with a concurrent reduction of the 
lattice distortion. 
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Fig. 5. Evolution of the transverse magnetization mo, the stag- 
gered magnetization ttIs, the uniform orbital order (/^) and 
the staggered orbital order parameter ts as a function of the 
applied magnetic field for three different sets of parameters 
ikBT/J,JG/K^) given by (i) (0.15,0.55), (ii) (0.15,0.45) and 
(Hi) (0.25,0.45). The points (i) — (iii) are also shown in the 
phase diagrams Figs. |3]and|31 



4.2 Applied magnetic field 

Now we introduce in our analysis the effect of an exter- 
nal magnetic field applied in the x-direction, that enters 
in the Hamiltonian by the term Timag given in Eq. 
For g^isH = 0.25 J we obtain the phase diagram shown 
in Fig. m Comparing Fig. [3] and [4] we see that the main 
effect of the application of the magnetic field is the dis- 
placement of the first order structural transition towards 
smaller values of G and, consequently, the reduction of 
the region of the phase diagram with AFM order, as ex- 
pected. Therefore, the lattice effectively becomes harder 
in the presence of a finite magnetic field. The first order 



nature of the structural transition is now present for a 
larger range of temperatures, up to the tricritical point 
TP shown in Fig. [D In addition, a critical end-point CE 
is now defined when the second order ksTAFM line meets 
the first order structural transition line. 



4.3 Metamagnetic transition 

The results of our calculations include a first order meta- 
magnetic transition (MMT). The characteristics of this 
transition are summarized in Fig. [51 Here we show, from 
panels a) to d), how the magnetic and orbital order pa- 
rameters change as a function of the field strength, for 
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different values of temperature T and elasticity G. The 
different (T, G) points are labeled by (?) — (Hi) and are 
also indicated in the phase diagrams Figs. [3] and HI The 
curves (ii) and (Hi) of Fig. [SK) show the discontinuous 
evolution of the magnetization towards a strongly polar- 
ized magnetic state by the application of a magnetic field. 
This magnetic transition is accompanied in our system 
by a structural transition where at the same metamag- 
netic critical field He an orthorhombic FO phase changes 
discontinuously towards an AFO phase with tetragonal 
symmetry, as shown in Fig. [5l;)-d). For the MMT to be 
observed it is not stringent that the zero-field phase has 
antiferromagnetic order, as long as it is close enough to 
such a phase. Notice however, that a MMT is only pos- 
sible for a soft lattice (small G) at low temperatures, as 
it is the case for the (T, G) points (m) and (Hi) of Fig. [3l 
For larger values of the elasticity, such as for (i), we find a 
continuos evolution of the order parameters. In summary, 
for the critical field Hc{T, G), we find the general behav- 
ior that both rising T or lowering G increases the critical 
field. 

Further insights concerning the conditions for the MMT 
to occur can be obtained from the iJ-G-phase diagram at 
zero temperature shown in Fig. [6l It is worth noting that a 
first order MMT (accompanied by a structural transition) 
can only occur at T = if we apply a magnetic field in the 
FO & AFM zone, and only for elasticity values belonging 
to the region G< < G < G>. For Eq < K/J, G< and G> 
are given by 



G< = 
G> - 



J{A + 4G) - AKeq ' 



4J {A [l/4:~27]l]~2B + C - 



(18) 




Notice that at G> there is a first order structural phase 
transition for T,H = 0. In Fig. [5] it can be seen that the 
metamagnetic critical field decreases as G is strengthened. 



Fig. 7. The orthorhombic distortion order eo — e= — ■^(/^) as 
function of the temperature for different values of the magnetic 
field between H = 0.2 and H — 0.5 (in units of J/gfis)- 



This qualitative behavior remains valid at finite but low 
temperatures. 



4.4 Magnetostriction and thermal expansion 



Since there is a close relation between structural and meta- 
magnetic transition the temperature dependence of the 
lattice parameters show a qualitative different low-temperature 
behavior for magnetic fields below and above the metam- 
agnetic transition. In our model, changes of the lattice 
parameters are considered by the orthorhombic distortion 
field e, related to the ferro-orbital order {I^) by Eq. p7|) . 
An increase of the c-axis is assumed to be proportional 
to Eq — s. Therefore, we show in Fig. [7] the tempera- 
ture dependence of Eq — e for different magnetic fields at 
G = 0A5K^/J. The T = critical field for this elasticity 
corresponds to He ~ 0.22 J/g/i^ as it can be deduced, for 
example, from the zero temperature G — H phase diagram 
of Fig. [6l For fields lower than He, a metamagnetic tran- 
sition will never be reached and the heating of the system 
by increasing the temperature drives the lattice towards a 
disordered FO phase through the crossover region (diffuse 
line in Fig. [3] and Fig. [l]). 

However, for fields H > He the system is in the meta- 
magnetic region. There is a low temperature AFO order 
and consequently, the distortion e is small. By heating 
the system, metamagnetism is destroyed and at the same 
time the lattice undergoes a first order structural transi- 
tion. This is seen in the e vs. T plot (Fig. [7]) by a jump 
of the distortion to a large negative value of Eq — e. If we 
keep heating the system, we reach again the FO region by 
passing the crossover zone. Obviously, the temperature of 
the l^*-order AFO/FO transition is larger for higher mag- 
netic fields. This is rcfiected in the evolution of the Eq — e 
discontinuity in Fig. [7]from H = Q.2bJ/giiB to 0.5J/gfiB- 
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5 Comparison to experiments 

Eventually, we will motivate our microscopic model in 
view of the experimental results. In particular, we will con- 
sider the phase diagram at zero-field as well as at the mag- 
netostriction and magnetization measurements that char- 
acterize the metamagnetic transition in Ca2-a;Srj:Ru04. 

5.1 Phase diagram 

The theoretical zero-field phase diagram (Figl3|) can be 
compared to the experimental T-x-phase diagram obtained 
in Ref. [5] and schematized in Fig.[TJ We relate the phases 
on the left-hand side of Fig. [3] (soft lattice) to the phases 
that are observed in region II of the experimental phase 
diagram. In both cases the symmetry of the lattice is re- 
duced. The tilting of the RuOe-octahedra observed in the 
real material is modeled by a distortion of the 2D lattice 
where the orbital and magnetic modes live. This reduction 
of symmetry is characterized by the strain field e or the 
ferro-orbital order (/^). For larger G (I^) is reduced and 
we find a transition to an AFO (predominantly) tetrag- 
onal phase where, in addition, FM order sets in at small 
temperatures. This is in fact the characteristics found in 
region III of the experimental phase diagram near x = 0.5. 
Therefore, it is not unfounded to relate the elasticity of the 
lattice (as modeled by G) to the Sr concentration x: both 
control the amount of distortion in their respective phase 
diagrams. Note that a Curie-like behavior of the orbital 
order (J^) « Ke/AksT is shown in Fig.[7l which accounts 
for the temperature dependence of experimentally mea- 
sured anisotropy of the spin susceptibility in the distorted 
(orthorhombic) region [5|25j . 

5.2 Metamagnetic transition 

The metamagnetic transition shown in Fig. [5] can be re- 
lated to the experimental MMT |2l3j in the following way: 
in zero magnetic field, at a temperature and doping (G 
in our language) leading us into the low-temperature FO 
zone of the phase diagram (region II of Fig[T]), we find a 
large strain e ^ eg and a zero-component of the magneti- 
zation along the x-direction. If we now turn on the trans- 
verse magnetic field, a finite component of toq appears, al- 
though for small enough fields the strain is still present in 
the system. For some critical field Hc{T, G) we find a first 
order transition in the magnetization mp which jumps dis- 
continuously to some larger value, while the strain drops 
simultaneously to e ^ Eq. The transcription of this to the 
experiments is that the octahedra returns to the structure 
it had before tilting. Also the c-axes adapts to the initial 
direction it had in region III of the phase diagram. The 
tetragonal symmetry needs, however, not to be restored. 
This behavior explains now the reversal of the structural 
distortion which occurs upon cooling at zero field, since 
applying a high magnetic field at low enough tempera- 
tures leads back to the old structure, as shown in Fig. [7] 
and seen experimentally [5]. 



This first order transition in the magnetization corre- 
sponds to the MMT observed in the experiments. Note 
that inhomogeneity is ignored in our description. What 
occurs as a discontinous first order transition here, would 
be a smooth crossover (MMT) when disorder, e.g. in al- 
loying Ca and Sr, is included [25] . 

The critical elastic constants G< and G> found in 
Fig. [HI bounding the segment on the G-axis where the 
first order magnetic and structural transitions are possi- 
ble at r = 0, can be mapped to the concentration values 
X of the experimental phase diagram that define the re- 
gion where the MMT can be observed. Therefore we may 
identify G< with x = 0.2 and G> with x = 0.5. This rela- 
tion is consistent with the experimental results that show 
a smaller energy scale for the transition at a; = 0.5 com- 
pare to the one at x — 0.2. The MMT is shifted towards 
lower fields when the Sr content grows from x — 0.2 to 
0.5. On the other hand, the first order structural transi- 
tion driven at zero temperature for G> in the theoreti- 
cal model can be related to the structure quantum phase 
transition of Ca2_a;Sr2;Ru04 at a; = 0.5. At G>, the meta- 
magnetic transition may be considered as occurring at 
zero-temperature and zero magnetic field, in agreement 
with the experimental measurements that shows that the 
MMT seems to be shifted towards a field close to zero for 
X = 0.5. 



5.3 Magnetostriction and thermal expansion 

The dependence of the first order magnetic transition on 
the temperature shown in Fig. [5] for points (ii) and [Hi) 
of the phase diagram, and the temperature dependence of 
the distortion of Fig. [7] can be compared to the experi- 
mental measurements, too. This is actually the expected 
behavior and reproduces some of the results shown in 
Ref. |2i3j . If we look for example at the magnetostriction 
measurements of Ref. [2|3I , where they show AL{H) / Lq 
along the c axis as a function of the applied magnetic 
field, these results can be interpreted as the response of 
the lattice structure to the metamagnetic transition. The 
jump of the magnetization is coupled to an increase of 
the c-lattice constant Loi since the structure of the lattice 
before the octahedra tilting is restored (see Fig.[5l:)). 

In addition, the thermal expansion coefficients and the 
integrated length changes measured below and above the 
critical field show that the pronounced shrinking of the oc- 
tahedra along the c-direction in zero field is successively 
suppressed by the field and turns into a low-temperature 
elongation at fields larger than He |2|3j . The experimental 
temperature dependence of AL/Lq for various magnetic 
fields, as a measure of the lattice distortion, follows a tem- 
perature dependence of e similar to the case of fields above 
and below He and shown in Fig. [T] In fact, for fields be- 
low the MMT the lattice evolves from the low-temperature 
FO order to the high-temperature PO region. On the other 
hand, for fields above the MMT the system shows a tetrag- 
onal AFO symmetry within the MMT- region which is 
suppressed by increasing the temperature, switching to 
an orthorhombic FO order. For higher temperatures the 



8 



R. Roldan, A. Riiegg and M. Sigrist: Interplay of metamagnetic and structural transitions in Ca2-a;Sra:Ru04. 



system looses orbital order and reaches the disordered PO 
region. 

6 Conclusion 

In summary, we have analyzed a microscopic model for 
the description of magnetic and structural properties in 
Ca2_xSra;Ru04which is based on the assumption that two 
of the three electron bands known in Sr2Ru04 are Mott 
localized in regions II and Ca-rich zone of region III (near 
X — 0.5) of the phase diagram |8]. The mean-field treat- 
ment of this model reproduces the basic magnetic and 
structural properties, as well as the metamagnetic tran- 
sition of this material. The elastic properties have been 
introduced assuming that the elastic constant depends on 
the Ca-concentration. In this way we draw the connection 
between our phase diagrams based on model parameters 
and the physical phase diagram and find a good qualita- 
tive and in parts quantitative agreement. Our most im- 
portant result is the magnetostriction effect in connection 
with the metamagnetic transition which agrees well on a 
qualitative level with recent experimental findings. 
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